Spiral Evolution in a Confined Geometry 
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Supported nanoscale lead crystallites with a step emerging from a non-centered screw dislocation 
on the circular top facet were prepared by rapid cooling from just above the melting temperature. 
STM observations of the top facet show a nonuniform rotation rate and shape of the spiral step as 
the crystallite relaxes. These features can be accurately modeled using curvature driven dynamics, 
as in classical models of spiral growth, with boundary conditions fixing the dislocation core and 
regions of the step lying along the outer facet edge. 

PACS numbers: 05.20.-y, 61.72.Ff, 68.35.-p, 68.37.Ef 



The control of mass flow in the fabrication and 
evolution of nanoscale structures is a key problem in 
nanoscience. Scanned probe techniques can be used to 
directly observe dynamical evolution, permitting the de- 
velopment of detailed models that can yield a predictive 
understanding of nanostructure evolution 
Here we address a classical problem — the facilitation of 
structural evolution by the presence of screw dislocations 
— in a novel context where the spiral step is confined to 
the top facet of a supported crystallite island prepared 
in a non-equilibrium configuration. 

Frank predicted that screw dislocations serve as a con- 
tinuous source of steps 0, crucial for step-flow growth 
on nominally flat terraces. In the presence of a super- 
saturation in the vapor, an infinitely long straight step 
starting at a dislocation rapidly winds into a spiral step. 
Once the curvature of the step at the dislocation core at- 
tains a particular critical value corresponding to the size 
of the critical nucleus in classical nucleation theory, the 
entire spiral rotates uniformly, with an angular velocity 
determined by the supersaturation of the vapor and the 
material parameters. Such growth has been studied by 
experiments H llol|. t heorv H [H E [3 and Monte 
Carlo simulations These studies of spiral evolu- 

tion have dealt with processes on infinite substrates. We 
show here how the dynamics of spiral steps in confined 
geometries can be appropriately modeled by accounting 
for the boundary conditions imposed by the edges of the 
nanocrystalline structures. 

We prepared three dimensional crystallites of Pb on 
a Ru(OOOl) substrate using methods described in previ- 
ous work 6, 15j. Direct observation of the relaxation of 
the crystallites following a rapid change in temperature 
was achieved using variable temperature scanning tun- 
neling microscopy (VT-STM) 0,0. By controlling the 
rate of cooling from 550 K to the measurement tempera- 
ture of 390 K, it is possible to create crystal populations 
where about 30% of the structures have a screw disloca- 
tion on the top facet that is embedded many layers into 



the crystal. Crystallites with such dislocations have been 
shown to relax to a final state much closer to the expected 
equilibrium crystal shapes (ECS) than those without dis- 
locations, suggesting that dislocations provide an easier 
kinetic pathway to crystal relaxation [ig. ITtI Ii9| . 

As shown in Fig. the shape of the step emerging 
from the dislocation qualitatively resembles the initial 
shape of a normal growth spiral. However, the later turns 
of the spiral merge into the edge of the crystallite instead 
of forming the expanding shape of a growth spiral. A ge- 
ometrical model describing the structure of the spiral is 
shown in Fig lb. Time resolved STM imaging of the 
structure as illustrated by Fig|2t shows that the spiral 
turns rapidly in the uphill direction, opposite to that of 
a growth spiral, corresponding to mass leaving the step 
edge. The time resolved measurements show strikingly 
that the shapes and angular velocities change substan- 
tially during each turn of the spiral step as successive 
layers of the crystallite peel off. In particular the motion 
of the step becomes so fast that only a few frames are cap- 
tured l20l| as the angle ^o, as defined in Fig^^, changes 
from TT to 27r. The radius of the top facet increases only 
very slowly compared to the motion of the spiral step 
and remains essentially constant for many turns of the 
spiral. To model this behavior, we assume, as in previ- 
ous models of spiral growth 0, 0| , that the step starts 
at a fixed point that corresponds physically to the core 
of the dislocation and moves according to the usual laws 
of motion for isolated steps. Here however the extent of 
the step is limited by the finite size of the island. We 
model this as a boundary condition confining the step to 
a fixed outer circle representing the island boundary. We 
assume there is no transport of atoms through the bulk 
and neglect the anisotropy of the step stiffness. The ex- 
perimental conditions and the low vapor pressure of Pb 
justify these assumptions. 

The local chemical potential fi of the curved step rel- 
ative to a straight step can be written down using the 
Gibbs-Thompson relation /j, = fls^K, where k is the 



FIG. 1: color on-line(a) STM picture of a 3-D crystallite of 
Pb on Ru with a screw dislocation leading to a spiral step on 
the top facet, (b) The top facet of a crystallite containing an 
embedded screw dislocation is modeled as a circle of radius 
R. The dislocation step is shown as the bold curve that be- 
gins at point P, which is a distance A from the center of the 
circle 0, and joins the circle smoothly at the joining point T. 
Beyond T, the dislocation step follows along the circle in the 
counterclockwise direction. The shape of the step for a given 
initial angle 9o is specified by the distance r as a function 
of the angle along the curve. The directions of the radial 
velocity of the step Vr and its component normal to the step 
Vn are separated by an angle 7. 



local curvature of the step, Qs is the atomic area and /3 
the step stiffness. The step velocity normal to the local 
direction of the step is denoted by w„ and can be de- 
rived from the chemical potential for different limiting 
cases. We use a model where the rate limiting process is 
attachment and detachment of atoms from the step re- 
gion ■ Then, the local step segment moves according 
to the chemical potential difference between the step and 
a reservoir chemical potential maintained by fast diffu- 
sion on the terrace. The local step velocity is given by 
Vii{k) = (Fyi/JlsfcT) [^j,{k) — /ires] whcrc Ta is the mobil- 
ity of the step edge and fires is the chemical potential of 
the reservoir. Substituting for Ta in terms of the rate 



FIG. 2: color on-line(a)STM snapshots of the time evolution 
of the step due to the dislocation, (b) corresponding snapshots 
at time t of the calculated evolution superimposed onto a 
single figure. Curve 1 (red) t — 0, curve 2 (green) t = (9/23)r, 
curve 3 (violet) t = (17/23)r, curve 4 (blue) t = (21/23)r, 
where T is the period of rotation. 



constant for detachment from the step edge kAD |22ll24 
the equilibrium concentration of adatoms on the terrace 
Ceq and the lattice constant a yields the following expres- 
sion for the local normal step velocity in terms of the 
local curvature: 



a'^CeqPkAD 

kT 



(1) 



In the above equation, we have expressed the reservoir 
chemical potential in terms of a cutoff curvature Kc- The 
atoms detaching from the dislocation step diffuse across 
the top facet to the edge of the crystallite, and down its 
sides, causing motion of the crystallite edge |2^. How- 
ever, the time scale of this motion is much larger than 
that of the motion of the dislocation step during each 
turn. Thus, to describe motion of the dislocation step on 
the top facet, the edge of the crystallite can be treated 
as an effective reservoir with Kc = 1/R where R is the 
average radius of the top facet. We assume that the spi- 
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ral step smoothly joins this outer circle, which with our 
choice of Kc does not move as the evolution proceeds. 

Figure n describes the coordinate system that we will 
use for our analysis. In these coordinates, the curvature 
is given by 



2r" 



(2) 



where the prime denotes a derivative with respect to the 
angle 9. The limiting angle 6*0 as the step emerges from 
the core has ^(^o) — and Eq. (0) then gives r'(6'o) = 
2/k(6'o). We can take r{0) = r'{e) ^ for < 6* < 6*0. 
Thus 6*0 is the smallest angle at which r'{9) ^ 0. Express- 
ing the normal velocity w„ = — (dr/dt) cos j in terms of 
the radial velocity yields the fundamental evolution equa- 
tion 



dr 
dt 



r2 + 2r'^ - rr" 1 
kT r [ (^2 + ^/2^3/2 - 

(3) 

This partial differential equation can be solved numeri- 
cally given an initial configuration. In order to reproduce 
the experimental results we assume the dislocation core 
is located a distance A = R/2 from the center of the cir- 
cle. We take a special initial condition where the step 
starts at the core with curvature Kc, extends to the circle 
with some arbitrary shape that joins the circle smoothly, 
and then follows the circle exactly for many turns in the 
counterclockwise direction as represented in Fig^ This 
choice of initial conditions fixes the both the dislocation 
core and the outer circle throughout the entire evolution, 
since there are no interactions between different step seg- 
ments in this model. The boundary condition confining 
the top step to the outer circle builds in the main ef- 
fect of step repulsions, which is the physical origin of the 
confinment. 

The numerical solution of the partial differential equa- 
tion yields a time dependent shape of the spiral that is 
shown in Fig.j^b. We find that after a brief transient pe- 
riod, the time-dependent shapes become independent of 
the initial conditions, and quantitatively reproduce the 
experimental observations. The marked change in evolu- 
tion of the step shape for between tt and 2t: is repro- 
duced by these solutions. In this range, the step shape 
contains highly curved regions as the step passes through 
the constricted geometry between the dislocation and the 
facet edge, leading to very large normal velocities in some 
parts of the step. 

The time-period of one rotation for a top facet of ra- 
dius 385 nm at a temperature of 390 K is observed to 
be 23 minutes. Substituting these values along with the 
value of the lattice parameter (a — 0.287 nm) into the 
model gives CeqkAoP/kT = 3.9 x 10^nm~^s^^ Using 
(3 = 0.315eV/nm at 390 K d HI we get c^qkAo = 



4.19 X 10^nm~^s~^. This is in good agreement with cal- 
culations from previous experiments |l8| . The fact that 
the numerical values are reasonable, along with the excel- 
lent agreement with the shape evolution, suggest that the 
additional mass transport mechanism, step edge diffusion 
(SED) does not play a major role in the overall evolution. 
This agrees with our physical expectation that for mass 
transport over large length scales, SED, as a "one-lane 
highway" , will be ineffective compared to mechanisms 
(AD or TD) in which material transport occurs across 
terraces 27J. 

In order to explore the connection with the growth 
spirals of Burton, Cabrera and Frank (BCF) Q, we con- 
struct approximate solutions called "uniformly rotating 
shapes" (URS). These shapes have, at any given instant 
corresponding to a value of j a constant angular veloc- 
ity (w = dO/dt) along the step from the core to the point 
T where they join the circle smoothly with a continuous 
first derivative r' . But unlike the case of a centered spi- 
ral, the angular velocity and the total arc length from 
the core to the joining point will change for solutions at 
different time instants. 

For such URS the left side of Eq. 13 equals r'u) and 
the equation reduces to an ordinary differential equa- 
tion (ODE). Its solution for r{9) for a give n angle 6q 
can be obtained by a shooting method |23| described 
here. We start with the known values of r(0o) = and 
r'{9o) = 2/kc, pick a particular value of r"(0) and nu- 
merically integrate the ODE. The value of r"(0o) is then 
adjusted until the solution joins the outer circle with a 
common tangent. The result for r{d) with the small- 
est absolute value of r" gives the URS for the given 
angle 6q. Since the angular velocity co is the same 
at different points along the step, we can calculate it 
from a Taylor expansion around r — with the result 
Lo = 3r"{eo)a*Ceq^kAD/SR^kT. FigureHshows the URS 
for the values of 6*0 in the time-dependent shapes. The 
URS joining the circle in the half plane above the diag- 
onal line in Fig. O closely resemble the time dependent 
shapes, and the times predicted using the instantaneous 
values of lu are consistent with the time dependent solu- 
tion. It can be mathematically shown that all the URS 
have to join the circle in this half plane. Here the distance 
of step segment from the core increases with increasing 
arc length along the curve, and a smooth URS satisfying 
the ODE can be found. 

However no such URS exist for a step joining the cir- 
cle in the half plane below the diagonal line, where there 
are regions of the step where this radial distance must 
decrease and r' is negative. If we try to define a lo- 
cal angular velocity for the time dependent solutions at 
each point along the step as the right hand side of Eq. O 
divided by r', then this estimate for the angular veloc- 
ity will diverge in regions when r' tends to zero and is 
negative when r' is negative. This shows that the time- 
dependent motion of the step in general cannot be well 
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FIG. 3: color on-line. Snapshots of the time dependent evolu- 
tion (bold red curves) and the URS at the same 60 (dashed 
blue curves) superimposed on each other, (a) The shapes cal- 
culated at min, 8 min, and 22 min closely match the URS 
for the same 60 whereas (b) those at 12 min, 18 min, and 
21 min do not resemble the URS at the same 60 ■ All URS 
have joining points lying above the diagonal dashed line that 
passes through the dislocation core and the circle center. 



approximated by purely rotational motion. 

The ODE is mathematically identical to the equation 
for spiral growth used by BCF. However, the physical 
interpretation of the cutoff radius is different from BCF. 
Unlike the BCF growth spirals, the shape and angular 
velocity of the URS are determined by the boundary con- 
dition at the crystallite edge. Thus the URS represents 
a different branch of the solution to the ODE. These so- 
lutions have a value of curvature at the core that is a 
local minimum, as opposed a global maximum as is the 
case for the BCF spirals. Unlike the expanding growth 
spirals, these solutions, when followed beyond the outer 
circle boundary, converge to a circle of radius R centered 
at the core. 

In conclusion, STM observations of mass transport via 
spiral step evolution can be effectively modeled using sim- 
ple geometric boundary conditions and the standard for- 
mulations of the continuum step model. The resulting 
motion is significantly different from previous theoretical 
and experimental studies of unconstrained spiral motion, 
revealing nonuniform angular velocity of the evolving spi- 
ral step arising from the asymmetry in the position of 
the dislocation core with respect to the confining geom- 
etry. Mechanisms of mass transport and structural rear- 
rangement in the confined conditions to be expected in 
nanoscale structures are likely to reveal many more such 
cases of interesting boundary effects. 
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